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ABSTRACT 

The vortices in a high-Tc superconductor with strong correlated pinning cen- 
ters have been studied numerically using the mapping to charged bosons in two- 
dimensions(2D) and the Monte-Carlo algorithm. Considering the viscous dissipa- 
tion of moving vortices we derived a nonlinear voltage response expression which 
describes different regimes and their crossover uniformly. This equation accords 
with experimental results. 
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I. INTRODUCTION 

The static and dynamic properties of vortices in the mixed state of high-Tc super- 
conductors(HTSC) have been intensively studied both experimently and theoretically in 
recent yearsGJ . For applying superconductors in external magnetic fields, it is important 
to minimize the dissipative loses of moving flux lines by improving the flux pinning with 
some kinds of strong correlated disorder (material inhomogeneities). The vortex dyns- 
mics with strong correlated pinning can be studied efficiently by exploiting the mapping 
between vortices and 2D bosons@. 

Similar to the physics of flux lines in a pure systemS, the statistical mechanics 
of vortices interacting with columnar pinning centers which are aligned parallel to the 
magnetic field may be mapped into the quantum mechanics of charged bosons in two- 
dimensions (2D)0. Table I summarizes the analogy between the vortices system, with 
the tilt modules e\ and thickness L(length of vortex), and the corresponding 2D charged 
bosons system®. 

Table I Boson analogy applied to vortex transport 
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In the Bose-Glass phase, the linear resistivity vanishes for low external current 
J <C J c , and the most important mechanics for vortex transport is "tunneling" between 
different columnar effect sites via the deformation of a pair of "superkinks" ^ . This 
is very closely related to variable-range-hopping(VRH) transport of charged carriers in 
disordered semiconductors!! , and leads to the highly nonlinear expression. 

By further exploring the analogy to two-dimensions(2D) Boses localized at randomly 
distributed defect sites, it leads to a " coulomb" gap in the distribution of the pinning 
energies g(e) near the chemical potential fi which separates the filled and empty energy 
levels®. In the limit of infinitely long-range interaction, A — ► oo, one would expect g(e) 
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to vanish near the chemical potential according to a power law 

g(e) = \e-p\ s . (1) 

This distribution affects the vortex transport properties. An inplane current J _L B 
induces a Lorentz force per unit length fi perpendicular to J, acting on all the flux 
lines: 

/ t = ^xJ (2) 

c 

which modifies the free energy of vortices system. In the boson picture this additional 
term represent an electric field E = z x J/c acting on the particles carrying charge <f)Q. 
In the spirit of the thermally assisted flux-flow(TAFF) model of vortex transport 01, the 
superconducting resistivity p = E/J may be written as 

pnp f exp[-U B (J)/kT] (3) 

where pf is a characteristic flux-flow resistivity, and Ub represents an effective barrier 
height which is of the type U B (J) = U (J /J) P M 

Besides this inverse power-law Ub(J), some other types have also been suggested, 
such as the Anderson-Kim model Ub(J) = U c (l — J c / J)&0 and the logarithmic barrier 

U B (J) = U c \n(Jo/J) ®- 

In a system with moving vortices, the total current density J = J s + J n , where the 
normal component 

J n = E{J)/p f (4) 



with pf = p n (B/B c2 ) which gives rise to dissipation and viscous dragtU. Thus the effect 
of Lorentz force on the effective barrier [/^(i.e., the saddel-point free-energy price for a 
flux line to leave its columnar pin) should be attributed to the supercurrent component 
of the total current density 

J s = J-E(J)/p f . (5) 
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This paper is orgnized as follows. In the subsquent section we briefly describe the 
model and the Monte-Carlo simulation procedure which have been succesfully used ®. 
In section III we derived the current-voltage characteristics with the consideration of 
the viscous disspation of moving vortices. In section IV we present results of simulation 
and compare them with some experimental results. Finally, a short summary concludes 
this work. 

II MODEL AND SIMULATION 

A.Model 

For describing a system which has No columnar defect sites randomly distributed on 
the xy plane, we use the model in Ref. Jj] with two-dimensional effective Hamiltonianfl 

j N D N D 

H = - niUjVirij) + ^ n A (6) 
and its grand-canonical counterpart 

No 

H = H-vJ2 n i- (7) 

i=l 

Here i,j = 1,2,..., No denote the defect sites. n« = 0, 1 represents the corresponding 

No 

site occupation number(nj = 1 if a flux line is bound in columnar defect rij), J2 n i is 

i=i 

the number of the flux lines. V(r) = 2eoKo(r / \) represents the repulsive interaction 
between the lines; the modified Bessel function Ko(r/X) describes a screened logarithmic 
interaction and Eq = (0o/4vrA) 2 with A the London penetration length. We have also 
included a random site energy t{, originating in the variation of pin diameters. Their 
distribution P can be chosen to be centered at < t >= 0, with width w. For simplicity 
we assume a flat distribution 

P(U) = 9{w - \U\)/2w (8) 

[Q(x) denotes the Heaviside step function]. 

For the interacting system®, we define single particle site energies Si as follows: 

dH Nd 

£i= «- = E n ^( r «)+*i- (9) 



For filled sites (n» = 1), Si is the energy required to remove the particle at sites i to 
infinity, for empty sites (rij = 0), correspondingly Si is the energy needed to introduce an 
additional particle from infinity to sites %. In the thermal equilibrium, the chemical po- 
tential \i separates the occupied and empty states. With the intervortex repulsion taken 
into account, the distribution of pinning energies g(s) can be viewed as an interacting 
single-particle density of states and may be obtained from the statistics of the energy 

levels £j. One would expect that the normalization of the energy distribution is 

+00 

/ g(e)de = N D /A = 1/d 2 (10) 
—00 

with A the area of the system. 
B.Simualtion 

Using a zero-temperature Monte-Carlo algorithm minimizing the total energy with 
respect to all possible one- vortex transfers^ , we reproduce the results of Ref . including 
the spatial configurations on the ground states and the distribution of the pinning energy 
g{s). We have performed extensive studies for the cases Np = 200, 400. In order to study 
the size effect, we simulate with Nd = 800 too. We have reproduced the result of Ref.Q 
of the size effect. Fig.l shows our result of the energy distribution, where one finds the 
" Coulomb" gap. 

Ill CURRENT- VOLTAGE CHARACTERISTICS 

Consider first the intermediate current regime J\ < J < J c (J\ = Uo/<fiod, and 
Uo =< Uk > the average of pinning energy 0) in which the motion of a single vortex 
is unaffected by the other vortex in the sampled. Driven by the external current J, 
a vortex will start to leave its columnar pin by detaching a segment of length Z into 
the defect-free region, thereby forming a half-loop of transverse size R. Considering the 
dissipation loses of moving flux lines, free-energy price of forming a half-loop of transverse 
size R by detaching a flux line segment Z into the defect-free region is apporixmately 

5F(R, Z) w e x R 2 /Z + U Z - f L RZ + E(J)<f> Q RZ/cp f (11) 
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From Eg. (JTT|) we estimate the saddle-point free-energe S(Fi)* and find the current- 
voltage relationship 

E{J)=p f exp[-(Ek/Kr)WJ.)] (12) 

where E k = d^eiUo and J s is the supercurrent component described in Eq.(|5|). 

For J L < J < J\(Jl is the current which satisfies R*(Jl) = L), we have to consider 
the configurational limitation imposed by other vortex. The most important thermally 
activated excitation will now be a double superkinkfi. This is the vortex analog of 
variable-range- hopping charge transport in disordered semiconductor @. The cost in free 
energy for such a configuration of transverse size R and extension Z along the magnetic- 
field direction will consist of three terms: (i)the double-superkink energy 2E k R/d stem- 
ming from the elastic term, (ii)the difference in pinning energies of the highest-energy 
occupied site, ~ p and the empty site at distance R with Si = p + A(i2), and (iii) the 
viscous dissipation of the vortex stemming from the motion of the current around the 
vortex kernel in the external mag netic field B, which is E(J)<f> Q RZ/p f EH. Thus the 
free-energy difference with respect to the situation without kinks and external current 



is 



5f w 2E k R/d + ZA(R) - f L RZ + E(J)<p Q RZ /cp f (13) 

— * 

The concentration available states as a function of R with D dimensions transverse to B 

At+A(K) 

(here D = 2) on the one hand equals dP j g(e)de, and on the other hand is simply 
given by w (d/R) D , thus A(i2) is to be determined from the equation@. 

J g(e)de = R~ D '. (14) 
ft 

Optimizing first for vanishing current J = gives the longitudinal extent Z* of the 
superkink function of its transverse size R*. 

2E k/d 



{dA/dR) R * 



(15) 
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Upon balancing the last term in Eq. ( |I3| ) against the optimized sum of the first two, one 
arrives at 

(J - E{J)/p f )<f>o/c = J s <P /c « A(R*)/R* (16) 

which through inversion yields a typical hopping range R*(J). Inserting back into 
Eq.(|I~3"D finally yields the result for the optimized free-energy barrier for jump, 



5F*(J) « {2E k /d)R*{J) (17) 
which we identify with the current-dependent activation energy in Eq.(^) 

E(J) « p f J exp[-(2E k /kTd)R*(J)]. (18) 
Considering a power-law form for the distribution of pinning energies 

g(e) = K\e-p\ s (19) 

one obtains 

5F*(J) = 2E k (J /J s y (20) 

where the transport exponent p 

P = 8 + 1 (21) 
v D+s+1 K J 

and the current scale 

Jo « c/<J)oK 1/{s+1) d 1/p . (22) 

Thus we get 

E(J) = Pf Jexp[-(E k /kT)(J /J s ) p ] (23) 

with J s described by Eq.(|^). 

Eq. (Hp can also be expressed in a general form@ 

y — x exp[— 7(1 + y — x) p ] (24) 
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with 




E(J) 

ULPf 

where Jhf = E{Jl)/ Pf which is much smaller than J^. 

Because of the size of the sample and the periodic boundry condition, the transverse 
size R*(J) cannot be greater than the sample width L. For J < J L , one would ex- 
pect that the current- vlotage characteristics will become Ohmic 0. This regime is also 
described by Eq. ([Zip . 

It is interesting to note, that the general form Eq. is also compatiable with other 
type of suggested Ug(J). For instance, an E(J) equation for type-II superconductors 
has been derived in connection to the Anderson-Kim model with replacing the total 
transporting current density J in Ub{J) = U c (l — J J J) with J s of Eq.(0)0, as 

0E(J) = Jp f exp[(-C/ c -W v + W L )/kT] (25) 

where W v = E(J) ■ B ■ A/pf is the viscous dissipation term of flux motion, Wl is the 
energy due to Lorentz driving force, Wl = J ■ B ■ A, the parameter A is a product of 
the volume of moving flux bundles and the range of force action. 
Eq.(|25|) can be expressed with a reduced form 

y = xexp[— 7(1 + y — x)] (26) 
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with 7 = U c /kT, x = Wl/U c and y = TV„/?7 C , which corresponds with Eq.(|2~3|) with 
p — 1. 

IV RESULTS AND DISCUSSION 

For simplicity, we assume that E^/k = 1. Thus temperature has unit E^/k. More- 
over, in this paper the energy is normalized by Eq, current is normalized by (pod/eoc and 
voltage is normalized by pf. Using the relations (|H|)-(|TE|) and the Monte-Carlo result 
g(s), we can numerically evaluate the current- volt age characteristics for any form of 
g(s). Results derived from the distribution of interacting pinning energies in Fig.l are 
shown in Fig. 2 which has highly nonlinear characteristics and a superconducting phase. 
The double-logarithimic plots of E(J) as the function of the current J in the inset of 



Fig. 2 is consistent with the result in Ref.[]I3 
From Eq. (|i8|) , we have 



a ^ ^> = -(2 W d jm. (27) 
a InJ aJ 



This immediately leads to 

~dR*(J) 



dR*{J) 



dJ 



{21 



Qmax dJ 

From Eq. (p!6|)-Eq. (p!8|) , we see that temperature will have little effect near the maxi- 
mal slope point. 

On the other hand, we get from the analytical equation Eq.( 



_ dlny 1 + pjx(l + y-xf 1 

D = — = 1 + cr = — r~T- ( ^iJ ) 

dlnx 1 + pjy(l + y — x) p 1 

At the maximal slope point (xj,^), we have 

p 2 j 2 xy(l +y- x) 2(p - 1] = l-p(x-y). (30) 

As we have mentioned in the previous section that the Eq. (JM|) is compatible with 
different types of suggested Ub{J), one would expect its general agreement with exper- 
imental results of various kinds of type-II superconductors. 

8 



In Fig. 3 we compare the current dependency of the slope a derived from Eq. 
and Eq.(pO|) with the wide range V ~ / data observed by Repaci et al. on the YBCO 
films!!!!. The agreement is rather well. 



We also compared our Eq.(24) with the scaling of isothermal E(J) curves observed 
by Koch et a/.0. Taking simple trial form 7 = 10(1 — T/T c ) s ■ T/T c with 5 = 0.5 
and p = 0.6, we get from Eq.(^) one hundred E(J) isotherms near the T g « 0.84T c (as 
observed in Ref.|lj]]) with T ranging from 0.74T C to 0.94T C . All the isotherms collapsed 
nicely onto two curves(T > T g and T < T g ), consistent with the scaling of v = 1.7, 



z = 4.8 as shown in Fig.4. The similar scaling result of Ref.|16| is shown in the inset of 
Fig. 4 which has the same scaling exponents of v — 1.7 and z = 4.8. 
V SUMMARY 

Using the Monte-Carlo simulation method introduced by Ref.[|J, we studied the 
current- volt age relation of type II superconductors which have columnar defect sites. 
Considering the viscous dissipation of moving vortices, we find that the three regime's 
current- volt age relation can be described by a unified equation Eq . ([24D . This equation is 
consistent with different types of Ub(J) so far suggested by different models and agrees 
with experimental observed E(J) data of wide range of temperature and current density. 
With proper trial effective barrier function Ub{T) it gives E(J) isotherms which can be 
collapsed onto two curves (T > T g and T < T g ) with scaling exponents found from 
experimental data. 
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Figure Caption 



FIG.l. Normalized distribution of pinning energies g{e) as function of the single-particle en- 
ergies e which is normalized by £q . Another parameters is w = 0.2. g(e) vanished near 
\x according to a power law \e — /i| 3 , thus the power p in Eq(19) is about 2/3. 

FIG. 2. The nonlineal current-voltage curve derived from FIG.l.-E(J). Voltage is normalized 
by pf while currrent is normalized by (ftod/eoc. 



FIG. 3. (a) The original experimental result of Ref.[15|. (b)The result cr/a max ~ ln(I/Ii) 
ploted from the results of Ref. fT|| are compared with analytical Eq. (f24|) and Eq. (|30|) (solid 
lines). The parameter p is taken as 1.8 and 7 ranges from 0.5 to 5.0 in 0.5 intervals. 
a{Ii) = a n 



' max ■ 



FIG.4. The collapsed data derived from Eq.(||), with v = 1.7, z = 4.8 and T g /T c = 0.84. 100 
curves are plotted for T > T g and T < T g . Inset is the original experimental result of 
Ref. 01 . 
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